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Abstract 

We present analytically the exact energy bound-states solutions of the 
Schrodinger equation in D-dimensions for a recently proposed modified 
Kratzer potential plus ring-shaped potential by means of the conventional 
Nikiforov-Uvarov method. We give a clear recipe of how to obtain an explicit 
solution to the wave functions in terms of orthogonal polynomials. The results 
obtained in this work are more general and true for any dimension which can 
be reduced to the standard forms in three-dimensions given by other works. 
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I. INTRODUCTION 



The important task of quantum mechanics is to find the exact bound-states solution of the 
Schrodinger equation for certain potentials of physical interest. Generally speaking, there 
are a few main traditional methods to study the exact solutions of quantum systems like 
the Coulomb, the harmonic oscillator [1,2], the pseudoharmonic [3,4] and the Kratzer [4,5] 
potentials. Additionally, in order to obtain the bound-states solutions of central potentials, 
one has to resort to numerical techniques or approximation schemes. For many of the 
quantum mechanical systems, most popular approximation methods such as shifted 1/N 
expansion [6], perturbation theory [7], path integral solution [8], algebraic methods with 
the SUSYquantum mechanics method and the idea of shape invariance, further closely with 
the factorization mrthod [9], exact quantization rule [10,11], the conventional Nikiforov 
and Uvarov (NU) method [12-25]. Some of these methods have drawbacks in applications. 
Although some other methods give simple relations for the eigenvalues, however, they lead 
to a very complicated relations for the eigenf unctions. 

The study of exact solutions of the Schrodinger equation for a class of non-central po- 
tentials with a vector potential and a non-central scalar potential is of considerable in- 
terest in quantum chemistry [26-35]. In recent years, numerous studies [36-40] have been 
made in analyzing the bound states of an electron in a Coulomb field with simultaneous 
presence of Aharanov-Bohm (AB) [41] field, and/or a magnetic Dirac monopole [42], and 
Aharanov-Bohm plus oscillator (ABO) systems. In most of these studies, the eigenvalues 
and eigenfunctions are obtained by means of seperation of variables in spherical or other 
orthogonal curvilinear coordinate systems. The path integral for particles moving in non- 
central potentials is evaluated to derive the energy spectrum of this system analytically [43]. 
In addition, the idea of SUSY and shape invariance is also used to obtain exact solutions of 
such non-central but seperable potentials [44]. Very recently, the conventional NU method 
has been used to give a clear recipe of how to obtain an explicit exact bound-states solutions 
for the energy eigenvalues and their corresponding wave functions in terms of orthogonal 
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polynomials for a class of non-central potentials [45]. 

Recently, Chen and Dong [46] found a new ring-shaped potential and obtained the exact 
solution of the Schrodinger equation for the Coulomb potential plus this new ring-shaped po- 
tential which has possible applications to ring-shaped organic molecules like cyclic polyenes 
and benzene. Very recently, Cheng and Dai [47], proposed a new potential consisting from 
the modified Kratzer's potential [48] plus the new proposed ring-shaped potential in [46]. 
They have presented the energy eigenvalues for this proposed exactly-solvable non-central 
potential in three dimensional (i.e., D = 3)-Schr6dinger equation through the NU method. 
The two quantum systems solved by Refs [46,47] are closely relevant to each other as they 
deal with a Coulombic field interaction except for a slight change in the angular momentum 
barrier acts as a repulsive core which is for any arbitrary angular momentum £ prevents col- 
lapse of the system in any dimensional space due to the slight perturbation to the original 
angular momentum barrier. 

The conventional Nikiforov-Uvarov (NU) method [12], which received much interest, has 
been introduced for solving Schrodinger equation [13-21], Klein-Gordon [22,23], Dirac [24] 
and Salpeter [25] equations. We will follow parallel solution to [47] and give a complete exact 
bound-states solutions and normalized wave functions of the D-dimensional Schrodinger 
equation with modified Kratzer plus ring-shaped potential, a Coulombic-like potential with 
an additional centrifugal potential barrier, for any arbitrary ^'-states using the conventional 
Nikiforov-Uvarov method. Our general solution reduces to the standard three-dimensions 
given by Ref. [47] in the limiting case of D = 3. 

This work is organized as follows: in section II, we shall briefly introduce the basic 
concepts of the NU method. Section III is mainly devoted to the exact solution of the 
Schrodinger equation in Z)-dimensions for this quantum system by means of the NU method. 
Finally, the relevant results are discussed in section IV. 
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II. BASIC CONCEPTS OF THE METHOD 



The NU method is based on reducing the second-order differential equation to a gener- 
alized equation of hypergeometric type [12]. In this sense, the Schrodinger equation, after 
employing an appropriate coordinate transformation s = s(r), transforms to the following 
form: 

+ ^n(s) + ^MS) = 0, (1) 

where a(s) and a(s) are polynomials, at most of second- degree, and f(s) is a first-degree 
polynomial. Using a wave function, ip n (s), of the simple ansatz: 

ipn(s) = (j) n (s)y n (s), (2) 

reduces (1) into an equation of a hypergeometric type 

oWJL*) + T(s)y' n (s) + Xy n (s) = 0, (3) 

where 

a W = *W§n> (4) 

(p'(s) 

r( s ) = f(s) + 2tt(s), t'(s) < 0, (5) 

and A is a parameter defined as 

A = A n = -nr'(s) - n{n ~^ ] a\s), n = 0, 1, 2, .... (6) 

The polynomial r(s) with the parameter s and prime factors show the differentials at first 
degree be negative. It is worthwhile to note that A or A n are obtained from a particular 
solution of the form y(s) = y n (s) which is a polynomial of degree n. Further, the other 
part y n (s) of the wave function (2) is the hypergeometric-type function whose polynomial 
solutions are given by Rodrigues relation 
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where B n is the normalization constant and the weight function p(s) must satisfy the con- 
dition [12] 

d t(s) 

—w(s) = —-w{s), w(s) = a(s)p(s). (8) 
as a { s ) 

The function ir and the parameter A are defined as 



\ 



'^mV_ 5(s)+Ms) , (9) 



A = A; + tt , (s). (10) 

In principle, since 7r(s) has to be a polynomial of degree at most one, the expression under 
the square root sign in (9) can be arranged to be the square of a polynomial of first degree 
[12]. This is possible only if its discriminant is zero. In this case, an equation for k is 
obtained. After solving this equation, the obtained values of k are substituted in (9). In 
addition, by comparing equations (6) and (10), we obtain the energy eigenvalues. 

III. EXACT SOLUTIONS OF THE QUANTUM SYSTEM WITH THE NU 

METHOD 

A. Seperating variables of the Schrodinger equation 

The modified Kratzer potential plus ring-shaped potential in spherical coordinates is 
defined as [47] 



Yir.fi) ■■■■ D, ( r -Hl) 2 + (11) 

r z sin v 



where f3 is positive real constant. The potential in (11) introduced by Cheng-Dai [47] 
reduces to the modified Kratzer potential in the limiting case of f3 — [48]. In fact the 
energy spectrum for this potential can be obtained directly by considering it as special case 
of the general non-central seperable potentials [45]. 



Our aim is to derive analytically the energy spectrum for a moving particle in the presence 
of a potential (11) in a very simple way. The D-dimensional space Schrodinger equation in 
spherical polar coordinates written for potential (11) takes the form [1,6] 



2/i 



1 a ^Ui 



-D-l 



dr 



dr 



]_d_ ( . d_\ _L d 2 2/^/3 cos 2 9 s 

/ 2 \ s^Ode \ 89 ) + sin 2 9 dp) 2 h 2 sin 2 9 



4>(r,9,<p) 



+ 



D P 



r — r. 



-E 



4>(r,9,<p) = 0, 



(12) 



where \i = ™^+^ 2 being the reduced mass of the two particles and ip(r, 9, ip) being the total 
wave function separated as follows 

^ £m (r, 9, ip) = R(r)Y e m (9, p), R(r) = r^^' 2 g(r), Y e m (9, p) = H{9)<S>{ip). (13) 



On substituting equation (13) into (12) leads to a set of second-order differential equations: 



1 d Ld-i± 
r ^ 1 dr \ dr 



T 2 



R(r) + | 



E-D P 



r — r f 



R{r) = 0, 



(14) 



l d ( . <£ 
Md9 V m d9 



m 



2/i(3 cos 2 9 



sin 2 9 h 2 sin 2 9 



+ £{£ + D-2) 



H{9) = 0, 



(15) 



d 2 §{p) 
d V 2 



+ m 2 $(p) = 0, 



(16) 



where L 2 D _ X = £(£ + D — 2). The solution in (16) is periodic and must satisfy the period 
boundary condition $(<p + 2n) = $(</?) from which we obtain 

1 



*m(V?) 



exp(±imp), m = 0,1,2, 



(17) 



Therefore, we are left to solve equations (14) and (15). After lengthy, but straightforward, 
calculations, Equation (14), representing the radial wave equation can be rewritten as [6]: 

d 2 g(r) 



dr 2 



+ 



^ (jp nu 4fiDere 1 V + Berlin 2 ) 



h 



h 2 



g(r) = 0, 



(18) 



where 



v = _(M - 1)(M - 3), M = D + 21. 



(19) 



The two particles in equation (18) interacting via Coulombic-like field have a slight change 
in the angular momentum barrier acts as a repulsive core which for any arbitrary I prevents 
collapse of the system in any space dimension due to the additional centrifugal potential 
barrier. On the other hand, equation (15) representing the angular wave equation takes the 
simple form 

m 2 + {2(1(3 /h 2 ) cos 2 Q 



d 2 H(9) cos 9 dH(6) 



dB 2 



sin# d9 



+ D-2) 



sin 2 9 



H{6) = 0. 



(20) 



Therefore, equations (18) and (20) have to be solved latter by using the NU method in the 
following subsections. 



B. The solutions of the angular equation 

In order to apply NU method, we introduce a new variable s = cos0, equation (20) is 
then rearranged as the universal associated-Legendre differential equation [45,49] 

d 2 H{s) 2s dH{s) v\\ - s 2 ) - m' 2 

" 1-s 2 ds + sin^l m = °' (21) 

where 

v' = £'(£' + D - 2) = £(£ + D - 2) + 2/2(3 /fi 2 and m' 2 = m 2 + 2/i(3/h 2 . (22) 

The solution of this equation has already been solved by the NU method in [45,47]. However, 
the aim in this subsection is to solve with different parameters resulting from the D-space- 
dimensions of Schrodinger equation. Upon letting D = 3, we can readily obtain the standard 
case given in [47]. Equation (21) is compared with (1) and the following identifications are 
obtained 

f(s) = -2s, a(s) = l-s 2 , a(s) = -v ' s 2 + v - m' 2 . (23) 
Inserting the above expressions into equation (9), one obtains the following function: 



7r(s) = < 



tt(s) = ± V /(z/' - A;)s 2 + A;-z/ + m /2 . (24) 

Following the method, the polynomial n(s) is found in the following possible values 

m's for k± — v' — m' 2 , 
—m's for ki — v' — m' 2 , 

(25) 

m! for ki — v'i 
—w! for k 2 = v' . 
Imposing the condition r'(s) < 0, for equation (5), one selects 

k\ — v' — m' 2 and tt(s) = —m's, (26) 

which yields form equation (5) 

r(s) = -2(1 + m')s. (27) 
Using equations (6) and (10), the following expressions for A are obtained, respectively, 

A = A n = 2n(l +rri) +n(n- 1), (28) 

A = i/-m'(l + m'). (29) 

We compare equations (28) and (29) and from the definition v' = £'(£' + D — 2), the new 
angular momentum £' values are obtained as 

£' = - (D ~ 2) + \\J{D - 2) 2 + 4(w + ym^+ 2/i/j/^ 2 )(n + 1 + ym 2 + 2/x/3/^ 2 ), (30) 

which can be reduced to the simple form 

£' = n + m', (31) 

in three-dimensions [47]. Using equations (2)-(4) and (7)-(8), the wave function can be 
written as, 

H m ,(0) = N em , sm(9) m 'Pi m '' m '\cos9), (32) 



where N^i m i = y ^ 2$+m'™ ^ * s ^ ne norma hzation constant given in [46,47] and 

n=-^^- + y(2e' + iy + 4t(D-3), 
with w! is defined by equation (22). 



(33) 



C. The solutions of the radial equation 



The solution of the SE for the modified central Kratzer's potential has already been 
solved by means of NU- method in [48] . Very recently, using the same method, the problem 
for the non-central potential in (11) has been solved in three dimensions by Cheng and Dai 
[47]. However, the aim of this subsection is to solve the problem with a different radial 
separation function g(r) in any arbitrary dimensions. We now study the bound-states (real) 
solution E < D e of equation (18). Letting 



/ 2/i 4fiD e r e „ 1 

e = \J-j2\ E - D e), " = —^l — , 7 = v+-ar e , 

and substituting these expressions in equation (18), one gets 

d 2 g(r) 



(34) 



-e 2 r 2 + ar — 7 



dr 2 



g(r) = 0. 



(35) 



To apply the conventional NU-method, equation (35) is compared with (1) and the following 
expressions are obtained 



r(r) = 0, cr(r) = r, a(r) = — e r + ar — 7 



Substituting the above expressions into equation (9) gives 



(36) 



7r(r) = - ± -y4e 2 r 2 + A(k - a)r + 47 + 1. 



(37) 



According to this conventional method, the expression in the square root be the square of a 
polynomial. Thus, the two roots k can be readily obtained as 



k = a ± £^47 + 1. 
In view of that, we arrive at the following four possible functions of 7r(r) : 

I + \er + + rl for h = a + eV 4 7 + 1, 



(38) 



7r(r) 



er + |v/ 4 7 + lj for fci = a + ey/¥f + 1, 
\ + er - |V 4 7 + 1 



(39) 



for k 2 = a — e^/47 + 1, 



er — \\[$rf + 1 for k 2 = a — ey/A^y + 1. 
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The correct value of n(r) is chosen such that the function r(r) given by equation (5) will 
have negative derivative [12]. So we can select the physical values to be 



k = a — E\J^ + 1 and ir(r) = - 



er--yj4<y + l 



(40) 



which yield 



r(r) = -2er + (1 + ^4 7 + 1). 



(41) 



Using equations (6) and (10), the following expressions for A are obtained, respectively, 



A = \ n = 2Ne, N = 0,1,2,..., 



(42) 



X = a-e(l + J47 + 1). 



(43) 



So we can obtain the energy eigenvalues as 



Em = D P 



2N + 1 + \J (M - 1)(M - 3) + 8(iD e r 2 Jh 2 + 



2 ' 



(44) 



where 



(M - 1)(M - 3) = W = (D - 2) 2 + 4£'(£' + D - 2) - 8^(3/h 2 - 1, 



(45) 



with £' defined in (30). Therefore, the final energy spectra in equation (44) take the following 
Coulombic-like form [7] 



E N i — D e — 



(N'T 



-, N' — 0, 1, 2, 



(46) 



where 



N' = 



1 



2N + \J {D — 2) 2 + 4f {£' + D - 2) + 8fi(D e r 2 e - (3)/h 2 + 1 



(47) 



is simply obtained by means of substituting equation (45) into (44). 

(i) If D — 3, equation (44), with the help of equation (45), is transformed into the 
following form 
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E!\[ n m D e 2 j 



8^Dyjn 2 _ 

^2N + 1 + \j(2n + l) 2 + 4m 2 + 4(2n + l)^m 2 + 2/2(3/ h 2 + 8/iD e rl/h 2 j 

(48) 

and it is consistent with equation (40) in [47]. 

(ii) If, D = 3 and (3 — (modified Kratzer potential), equation (44) is transformed into 
the form 

E n = D e ■ gggg ^ . (49) 

(l + 2n + v /l + 4£(£+l) + 8/i J D e r 2 /^ 2 ) 

and it is consistent with equation (14) in [48]. 

Let us now turn attention to find the radial wavefunctions for this potential. Using r(r), 
7r(r) and a(r) in equations (4) and (8), we find 

4>{r) = r (v^+T+i)/2 e -^ ; (50) 

p (r) = r Vmi e ^er^ (51) 

Then from equation (7), we obtain 

rfN 

y n (r) = B n r-^e 2 ^^ (r"+V^ e -^ ? (52) 

and the wave function g(r) can be written in the form of the generalized Laguerre polyno- 
mials as 

g(r)=C N , L r L+1 e-* r L 2 N L+1 (2er), (53) 

where L can be found easily from equation (47). Finally, the radial wave functions of the 
Schrodinger equation are obtained 

R(r) = C N , L r L ^ D -VI 2 e- £r Lf +l (2er), (54) 

where 
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e = 



(55) 



with N' is given in equation (47) and Cn,l is the normalization constant to be determined 



below. Using the normalization condition, / R 2 (r)r dr = 1, and the orthogonality re- 

o 

lation of the generalized Laguerre polynomials, / z n+1 e~ z [L^(z)] 2 dz = ( 2w + T /+- [ X ra + ? ?) ! ; we 
have 



C 



N,L 



(2e) 2L + 3 N\ 



\ 2(N + L + l)(N + 2L + iy: 
Therefore, we may express the normalized total wave functions as 



(56) 



(2e) 2L+3 (2£' + l)(f -m')\N\ , „ , 
- i M¥+ Ly.(N + L TWN + 2L + 1)! eM-cr)Lf+H2er) 



x sin(^) m 'Pi m '' m ')(cos^) exp(±imv?). 



(57) 



IV. RESULTS AND CONCLUSIONS 

In this paper, the Schrodinger equation in any arbitrary dimensions has been solved for 
its exact bound-states with a recently proposed modified Kratzer potential plus ring-shaped 
potential by means of a simple conventional NU method. The analytical expressions for the 
total energy levels of this system is found to be different from the results obtained for the 
modified Kratzer's potential in [48] and also more general than the one obtained recently in 
three-dimensions [47]. Therefore, the noncentral potentials treated in [45] can be introduced 
as perturbation to the modified Kratzer's potential by adjusting the strength of the coupling 
constant (3 in terms of D e , which is the coupling constant of the modified Kratzer's potential. 
In addition, the angular part, the radial part and then the total wave functions are also 
found. Thus, the Schrodinger equation with a new non-central but separable potential has 
also been studied (cf. [45] and the references therein). This method is very simple and useful 
in solving other complicated systems analytically without given a restiction conditions on 
the solution of some quantum systems as the case in the other models. Finally, we point out 
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that these exact results obtained for this new proposed form of the potential (11) may have 
some interesting applications in the study of different quantum mechanical systems, atomic 
and molecular physics. 
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